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Abstract 

A perturbative study of the Schrodinger equation in a strong elec- 
tromagnetic field with dipole approximation is accomplished in the 
Kramers-Henneberger frame. A prove that just odd harmonics ap- 
pear in the spectrum for a linear polarized laser field is given, assum- 
ing that the atomic radius is much lesser than the free-electron quiver 
motion amplitude. Within this approximation a perturbation series is 
obtained in the Keldysh parameter giving a description of multipho- 
ton processes in the tunneling regime. The theory is applied to the 
case of hydrogen-like atoms: The spectrum of higher order harmonics 
and the above-threshold ionization rate are derived. The ionization 
rate computed in this way determines the amplitudes of the harmon- 
ics. The wave function of the atom proves to be rigid with respect to 
the perturbation so that the effect of the laser field on the Coulomb 
potential in the computation of the probability amplitudes can be 
neglected as a first approximation: This approximation improves as 
the ratio between the amplitude of the quiver motion of the electron 
and the atom radius becomes larger. The semiclassical description 
currently adopted for harmonic generation is so rederived by solving 
perturbatively the Schrodinger equation. 



* e-mail: marcofrasca@mclink.it 



1 



Availability of powerful sources of laser light has permitted, in 
recent years, the realization of experiments through gaseous media 
that have shown several new physical effects as photoionization with 
a number of photons absorbed by the electron well above the ioniza- 
tion threshold and generation of a broad range of harmonics of the 
laser frequency Q . This latter effect could have a lot of technological 
applications and, as such, has been widely studied both theoretically 
and experimentally. 

The possibility to turn a physical effect into a practical applica- 
tion is strongly linked with the availability of a satisfactory theoretical 
model. But, it is common belief that, due to the intensity of the laser 
field, no perturbation theory can be done. The main aim of this pa- 
per is then to show how perturbation theory can be straightforwardly 
applied also for intense laser fields and analytical expressions can be 
computed for any kind of multiphoton process, at least for hydrogen- 
like atoms. The development parameter turns out to be the the square 
root of the ratio between the ionization energy Ib and the ponderomo- 
tive energy Up proportional to the intensity of the laser field, known in 
literature as the Keldysh parameter 7. The regime of a small Keldysh 
parameter characterize the so-called tunnelling regime that is the one 
of interest here. 

Theoretical approaches to multiphoton processes are non-perturbative 
in nature and resort to Floquet theory as in , numerical methods ap- 
plied directly to the Schrodinger equation as done firstly in |^ or semi- 
classical models [Q]. On the basis of the semiclassical ideas, a quantum 
theory for harmonic generation has been obtained by L'Huillier and 
coworkers in j^: Our theory permits to justify the main assump- 
tions of the quantum theory of these authors, so that, in turn, the 
semiclassical ideas prove to be a fairly good description of harmonic 
generation. 

The approach we apply to the Schrodinger equation for an atom 
in an electromagnetic field can be easily understood using a two-level 
model, widely used for harmonic generation Q. This model has the 
Hamiltonian (here and in the following we will take h = c = 1) 

H = -^a^ + ncos{ujt)ai (1) 

being ujq the level separation, 0, the intensity of the laser field and uj 
its frequency, ai and (T3 are Pauli matrices. If is small with respect 
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to ujQ, standard perturbation theory applies by interaction picture 
through an unitary transformation that removes the unperturbed part 
of the Hamiltonian: This gives a Dyson series in the smah development 
parameter 0,/{ujo ± w), out of resonance. Recently, duality has been 
introduced in perturbation theory and a dual interaction picture 
has been devised where one does an unitary transformation to remove 
the perturbation. For the above Hamiltonian one has to take U = 

e ^'^^ that yields the transformed hamiltonian ||8| 

Hp = ye "0 ^ Vs 

where now, perturbation theory can be done for Vl ^ loq^lo. We 
see straightforwardly that the unperturbed part of the Hamiltonian is 
"dressed" by the laser field and so, the energy levels are shifted. Then, 
the perturbation has odd and even harmonics of the laser frequency 
and both can appear in the spectrum. But, probability amplitudes 
that enters in the computation of the spectrum do not depend on 
the unitary transformations one does on the Hamiltonian and the 
states. So, we have sketched the physics of the two-level model in 
an intense monochromatic field just through dual interaction picture. 
Although, as we will show, the two- level model does not apply for 
current experiments with atomic samples as in this case one observes 
just odd harmonics in agreement with our full theory and it is not 
just a problem of a proper experimental setup, nevertheless it could 
have a wide range of applications in magnetic resonance experiments, 
for some other kind of media as optical cavities P| or wherever the 
conditions one meets for atomic samples are no more fulfilled. 

The dual interaction picture applies in the same way also to the 
Schrodinger equation in a semiclassical laser field and in the dipole 
approximation, as currently treated in literature The correspon- 
dence with the two-level model above is remarkable. The Hamiltonian 
in this case is 

H = ^ + y(x) + ^A(t)p + ^A\t). (3) 
Zm m 2m 

By the unitary transformation U{t) = exp {^-i^ Jq dt'A{t')p - dt'A^{t') 
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the above Hamiltonian transforms into 

Hkh = U\t) + y(x)^ U{t) = 1^ + y[x + a(t)] (4) 

being a(t) = Jq dt'A{t'). This is the well-known Kramers-Henneberger 
Hamiltonian and the unitary transformation above define the so called 



Kramers-Henneberger frame [10| that shows as the effect of the elec- 
tromagnetic field is to introduce a time-dependent translation on the 
potential of the unperturbed Hamiltonian by a length a(t). The laser 

E(t) 

field can be modeled as A(t) = , fxcosfcjt) -|- sin(u;t)l for a 

general ellipticity parameter ^. Here we consider the simplest case of 
a linear polarization ^ = and an instant rising of the laser field, that 
is E{t) = const. So, one has Q 

Hkh = ^+ — " "^'^ ^' + V i' [e-^-^+(-l)^e-'^->,(x) 



2m J_Ai TT . A2 

with 



(5) 



f^i- dx' (x^) 
Vk{x) = — V{x-x',y,z) ' = (6) 

being T\^{x) = cos(/i: arccos(3;)) the /c-th Chebyshev polynomial of first 

kind and Ar = = \l the maximum free-electron quiver mo- 
tion excursion. This length is pivotal in the study of atoms in an 
intense laser field as generally one has Al » a, being a = the 
Bohr radius. One can see that, as for the two-level model, we have 
the potential of the unperturbed part of the Hamiltonian "dressed" by 
the laser field and all the harmonics, odd and even, are present in the 
perturbation. We can now show that, in all the current experiments 
where the potential l^(x) depends just on r = |x| and 3> a, being a 
the Bohr radius of the atoms in the sample, then just odd harmonics 
appear in the spectrum. Indeed, we can rewrite eq.(^ as 

t-fc(x) = /' dx'Y{^{x-\j^x'f^y^-rz'') (7) 

If the laser field is enough intense, a series in ^ is obtained if one 
develops eq.(|^) in Taylor series as 

Vk{x) = p^dx' V{^{x-XLx'f + y^ + z^^ 



x=0,y=0,z=Q VrVT" 



X 
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-X 



\x'V\X,\x'\)^-^^l= + .... (8) 

-1 \x'\TTyJ\ - x'^ 



Despite its appearance, the terms of this series can be evaluated for a 
Coulomb potential and proved to be finite assuring the convergence. 
This is due to the fact that in this case the integrals can be computed 
analitically. Then, from the above expression two main conclusions 
can be drawn. Firstly, multiphoton effects are due to a dipole induced 
on the atom in the same direction as the electric field of the laser 
and secondly, Chebyshev polynomials have a definite parity and due 
to the symmetrical range of integration, only odd polynomials give 
a non-null contribution to the second term, while the first term has 
no physical consequences and in the following will be neglected. So, 
only odd harmonics contribute to the spectrum while, even harmonics 
are quadrupole radiation and then strongly depressed. Indeed, for a 
Coulomb potential one obtains 

r Ze"^ 
V2n+i{^) ~ -i{-ir — {2n + 1) — . (9) 

This result, that does not involve any other approximation beside the 
simmetry of the potential and the amplitude of the quiver motion 
of the electron with respect to the atomic radius, supports in some 
way the physical view recently given in |11], where it is assumed that 



the electron recolliding with the atomic core, emits bremsstrahlung 
radiation that is cut off at the maximum amplitude of the quiver 
motion of the electron, producing in this way just odd harmonics. 

To complete the above discussion before introducing perturbation 
theory, we have to study the "dressed" potential vq. This should be 
managed differently from the time-dependent part. Indeed, we have 
to separate the original potential V{r) from the shifts induced by the 
laser field on the energy levels of the atom. This can be obtained by 
a Taylor expansion as 



, dx' V(J(x - Xlx')'^ + + 



1 vr Vl - 



= V{r) + 5LV{-K) 

= V(r) + -^ V'(r)y^ + V'(r)z^ + V"(r)x^r + • • • (10) 
4r'^ L J 

where is seen that only even terms survive and higher order terms fall 
off very rapidly with r. The above expression assumes a very simple 
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form for a Coulomb potential 



t;o(x) = 



n=l 



(11) 



being An = ^\ dxx^^ / {tt^/ 1 — x"^) and P„ the n-th Legendre polyno- 
mial. This way to express the dressed Coulomb potential gives us a 
way to prove that the wave function is "rigid" with respect to the per- 
turbation using standard Rayleigh-Schrodinger perturbation scheme, 
for the kind of problems we discuss here. But, it should be pointed 
out that for stabilization things are quite different |l2|| . 
The equations for the amplitudes are given by 



la. 



+00 

i'=a„(t)(m|?;fc(x)|n)[e-^(^"-^'"-'='^)* + (-l)^e-^(-^"--^™+tl^^] 

n k=l 

having set En = En + {n\6LV{x)\n), being 6lV{x) the part of the 
static potential due to the laser field. At this point, all the machinery 
of standard perturbation theory applies [^]. For our aim, we have to 
show that the Rayleigh-Schrodinger part gives indeed a small contri- 
bution to the amplitudes. By assuming the atom initially in its ground 
state, this contribution is 



E^ — En 



-'m 



Using eq.(pT|) is easy to verify that no contribution comes for m = 2 
as (2|(5il/(x)|l) = but the degeneracy of level 2 is removed by 
the dressed potential as one has {m = 2,1 = 1,1^ = 0\5LV{x.)\m = 
2,1 = 1,1^ = 0) = ZeV(240a)(AL/o)2 and (m = 2,/ = = 
±l\6LV{x)\m = 2,1 = l,h = ±1) = -Ze^/{ma){\L/af while 
(m = 2,/ = 0,/^ = Q\5LV{-K)\m = 2,Z = 0,/^ = 0) = 0. Indeed, 
one can see that all the states having m even do not give a first order 
contribution even if the level shift is not null, while the level-shift is al- 
ways when I = 0. Instead, for m = 3 one has e.g. {m = 3,1 = 2,1^ = 
0\5LV{x.)\m = 1,1 = 0,1^ = 0) = Ze2^/T50/(10800a)(AL/a)2 and for 
the level shifts {m = 3,1 = 2,1^ = 0\6LV{x.)\m = 3,1 = 2,1^ = 0) = 
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ZeV(5670a)(AL/a)2 - Ze^/ {136080a) {Xi/a)^ and {1\6lV{x)\1) = 0, 
so the correction of eq.([l3|) turns out to be 



/ISO 



(^) 



4Ut) - — , r°V"\ .4 (^^^^^^ - 1) (14) 



4 I 

9 ^ 5760 



m \ a J 136080 V a / 



that is indeed neghgeable and the wave function turns out to be "rigid" 
with respect to the deformations introduced by the laser field. This 
is even more true as larger become the ratio A^^/a. The reason for 
this is that only a finite number of terms of eq.(0) give a non-null 
contribution to the matrix elements. It is interesting to note that for 
stabilization of an atom in intense laser field the situation is exactly 
the contrary as one should be able to diagonalize the Hamiltonian 
Hq = p'^/2m + fo(x) being the time-dependent part negligible, an 
approximation that becomes exact in the limit of infinite frequency of 
the laser field H, 0]. 



Then, the iterative procedure to solve eq.(12) can be applied to 
compute the probability transition for any process. This approach 
implies that off-resonant contributions should be sistematically ne- 
glected. In this way, a golden rule is straightforwardly obtained as 

+ 00 

Pi^f = 27rJ2 K^bn(x)|/)p5 [Ef - Ei - nu;] (15) 

n=l 

from which several results for multiphoton processes can be obtained. 
It is assumed a continuum of final states to sum over so that excited 
levels can decay, otherwise quantum resonance theory applies [|l5|] and 
Rabi hopping is obtained. In any case, going to second order gives 
a. C.Stark shifts of the energy levels. Rabi frequency due to resonance 
with the k-th harmonic of the perturbation with two levels m and n 
of the atom is (ref.ll|) ^ = |(m|ufc(x)|n)|. 

From eq.(|T5|) we can easily compute the rate of above threshold 
ionization. For hydrogen-like atoms [eq.(|^] and assuming the atom 
initially in its ground state one has 



32 cj^ 2 

TJ 



I 



B 



{2n + l)uj 



1 



{2n + l)uj 



(16) 



being no the minimun integer for which (2no + 1)uj — Ib > 0. It 
has been used the fact that, as shown above, for the ground state of 
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hydrogen-like atoms there is no shift by the part of the static potential 
due to the laser field, that is (l|(^Ly(x)|l) = for Coulomb potential. 
Beside, a plane wave is assumed for the particle in the final state 
to make computation simpler. By taking ref. |14| for experimental 
results , we can check the above expression for helium and neon that 
show a large plateau in the tunneling regime. So, we have Up = 155 
eV being the intensity 1.5 x 10^^ W/cm^, uj = 1.177 eV and Ib = 
24.59 eV. Then, 7 .4 and T « 0.026 eV, that is small as it should 
be expected. The same computation for neon gives approximatively 
0.02 eV. 

To analyse the question of harmonic generation, one has to com- 
pute < X >= {t)\x\^ {t)) . To complete this computation, we assume 
that no intermediate resonance is present and will justify this assump- 
tion a posteriori through the quantum resonance theory of ref.[^], 
that here applies. So, let us take an atom initially prepared in its 
ground state as to have ai(0) = (5jo- From eq.(p^) one has 

Eq — Em — «e 



+00 

Y,i^an{t){m\vk{-^)\Q) 

k=l 



-i{EQ~Em — kiO—it)t p—i{Eo~Em+ku)—ie)t 

+ (-1)'=-^ . ^17) 



Eq — Em — kiv — ie Eq — Em + kuj — ie 



with the limit e ^ understood as to have -^^^ = =fiTT6{x), being 
P the principal value. As is customary in perturbation theory, we 
keep just those terms that are near resonant with the harmonics of 
the perturbation: The only possibility left is the continuous spectrum, 
as it should be with the current understanding of harmonic generation. 
So, we take 

+00 i{Ep-Eo-kuj+ie)t 

«P W « - E i'(pb.(x)|0)(-l)^-^ ^ (18) 

being p the momentum of the particle in the continuous part of the 
spectrum. Now, we specialise this expression to the case of hydrogen- 
like atoms having 

Vg2 +00 i(^Ep-Eo-(2n+l)u+ie)t 

"pW - ^ gP" + l)(Pl^l°> E,-E„- (2n + 1). + u ■ 

Then, for the dipole moment one has 

< X E ap{t)e-'^^^-^°^'{0\x\p) + c.c. (20) 
p 
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After passing from the sum to integration through V J ^^^^3 and 

taking for the final state a plane wave, one gets the final expression 
for the harmonic spectrum 



3 

+00 ™2 



64 ZC iO C V A Jun / / ^ N / X 

< ^ >^ -i:i-f?r^ E 7 ;;T5 sm((2n + l)i^t) (21) 



32 Up n=no 



being x„ = ^^""^3(7^ ■ The normalization to 3{7p for originates 
from the fact that, from the above expression, the intensities of the 
harmonics reduce as the factor 3Up increases. Then, if the Keldysh 
parameter 7 is enough small we can take 

< ^ >^ E — + 1)^*) (22) 

32 p n=no Xn 

SO that, only for Xn < 1 the harmonic amplitudes are large. This is the 
approximate cut-off law found out through semiclassical methods in 
ref.Q. It should also be stressed the existence of a minimun harmonic 
order ng that should be expected due to the close connection between 
harmonic generation and multiphoton ionization. Indeed, this lower 
bound comes out from the phase space through the integration of the 
Dirac function both for the golden rule ( |l5|) and for the computation 
of the dipole moment < x >. Then, one gets = 10 and 9 for helium 
and neon respectively, that means harmonic 21 for the starting point 
of the spectrum in the regime of interest. It should be pointed out 
that the above equation for < x > has to take properly into account 
the ionization rate F of eq.(16) as to have at last 

+00 - 

< ^ >- -^^7^ E . "^"..5 «in((2n + X)u.t)e-^K (23) 



02 "-^p n=no 



[xn + Ij ) 



One can estimate the constant factor that determines the amplitude 
of the harmonics ^^rrr^l^- Indeed, for helium one obtains approxi- 

mately .32 x 10~^ eV~^ and for neon about .12 x 10"'' eV~^, showing, 
as it should be, a larger amplitude for neon. 

A further analysis concerns the effect of intermediate resonances 
on the spectrum of harmonics. On the basis of the theory of ref. [^], 
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one can write down eq.(|T^) as 



aUt)^-J2iHp\vkir)m-l) 



^i{Ep — Eo—kLij+ie)t 



( 



^t) (24) 



COS 



k=l 



Ep — Eq — kuj + ie 



with fi/j the Rabi frequency computed taking in account the reso- 
nances between the ground state and other discrete levels. To compute 
the above expression we assumed that the atom is initially prepared 
in its ground state so that, ao{t) = cos (J^tj, essentially the rotating 
wave approximation. It is easy to realize that one gets the harmonics 
in the spectrum shifted by the quantity 

The theory above could have wide applicability as, in principle for 
any multiphoton process one is able to compute analytical formulae 
to compare with experimental results. For instance, an improvement 
easy to implement is to use a full Coulomb wave function also for the 
final state in the above computations. On the other hand, even if 
major features of multiphoton processes are described by this theory, 
several problems are surely opened up as the applicability of the theory 
for an ellipticity parameter ^ 7^ 0, the introduction of a slower rising 
of the laser field or how to take into account all the features that real 
experiments have for harmonic generation. Beside, when the inten- 
sity of the laser field becomes too high the above approach should be 
properly modified as relativistic effects enter into the physical picture 
and, e.g. even harmonics can also be significant jl^ . Experiments to 
generate even harmonics are also carried out using solid surfaces as 
in [1^. Anyhow, it should be stressed how the possibility to derive a 
perturbative solution to the Schrodinger equation could give a chance 
to check models of multiphoton physics that no other approach offers. 
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